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\ 1. Introduction 

' Mathematical techniques could be developed into a more suitable and significant course by pre- 

senting them within the more general Sturm-Liouville theory in L2. The Sturm-Liouville Problems 
are important in many areas of science, engineering and mathematics. It is known that the spectral 
characteristics are spectra, spectral functions, scattering data, norming constants, etc. According 
to the theory a linear second-order differential operator which is self-adjoint has an orthogonal se- 
quence of eigenfunctions in L2. Spectral properties of Sturm-Liouville operators are often derived, 
directly or indirectly, as a consequence of an established link between large distance asymptotic 
behavior of solutions of the associated differential equation and spectral properties of the corre- 
sponding differential operator. Sturm-Liouville problems are divided into regular and singular type. 
\ Differential equations such as Bessel, hydrogen atom, Hermitte, Jakobi, and Legendre equations 

can be transformed into Sturm Liouville equations. There are many studies on these issues, [1 — 5]. 

Fractional calculus is "the theory of derivatives and integrals of any arbitrary real or complex 
order, which unify and generalize the notions of integer-order differentiation and n— fold integra- 
tion" [6 — 13] . In recent years, fractional calculus has been of great interest and there has been 
many results on the topic. It has been proved that many systems in different fields of science and 
engineering can be modeled more accurately using fractional derivatives [14 — 15]. We note that or- 
dinary derivatives in a traditional Sturm-Liouville problem are replaced with fractional derivatives 
and the resulting problems are solved using some numerical methods [16, 17]. Furthermore, Klimek 
and Agrawal [18] define a fractional Sturm-Liouville operator, introduce a regular FSLP and inves- 
tigate the properties of the eigenfunctions and the eigenvalues of the operator. In this paper, our 
purpose is to give regular fractional Sturm-Liouville problem for diffusion operator (FSLPDO) and 
prove fundamental properties of spectral data for the operator. 

2. Preliminaries 

Definition 1. Let Re (a) > 0. The left-sided and respectively right-sided Riemann-Liouville 
integrals of order a are given by the formulas 

X 

a 
b 

1 I > \a-l 



{Ib,-f) (x) = p^y / ( s " x ) f (s) ds, x<b (2) 



where T denotes the gamma function. 

Definition 2. Let Re (a) > 0. The left-sided and respectively right-sided Riemann-Liouville 
derivatives of order a are defined as 

(*) = D (l£°f) (x) x>a (3) 

(s) = -£> (/Jr/) (*) x<6 (4) 
Analogous formulas yield the left- and right-sided Caputo derivatives of order a: 

( C D« + f) (x) = (/£?£>/) (x) x>a (5) 

( c D%_f)(x) = (l 1 b Z a (-D)f)(x) x<b (6) 

Definition 3. [13] The general function p ^ q (z) is defined for z G C, complex ai,bj G C, and 
real ai,/3j 6R(1 = l,...,p;j = l,...q) by the series 
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_oo n r (°z + 

L (&i,/3i)i P 



( a i>«i)i,p . =J^-B £1 ( 7 ) 

fc=o ft r(6j + /3-A;) fc! 



,9 



where z, a;, G C, ai, /3j G R (/ = 1, ...,£>; j = 1, ...(?) . This general Wright function was investigated 
by Fox who presented its asymptotic expansion for large values of the argument z under the 
condition 

q p 

J>-I>>1. 

j=l 1=1 

If these conditions are satisfied, the series in (7) is convergent for any z G C. 
Teorem 4. [13] Let ai,bj G C, and au,/3j el(i = l, ■■■,P',j = 1, •••<?) and let 
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I) If A > —1, then the series in (7) is absolutely convergent for all z G C. 

II) If A > —1, then the series in (7) is absolutely convergent for |z| < 5 and for \z\ = <5 and 
3?(/i) > 1/2. 

Definition 5. [19] Let (X, d) be a non-empty complete metric space. Let T : X — >■ X be a 
contraction mapping on X, i.e.: there is a nonnegative real number q < 1 such that 

d(T(x),T(y))<qd{x,y) 

for all x, y in X. Then the map T admits one and only one fixed-point x* in X (this means 
T(x*) = x*). Furthermore, this fixed point can be found as follows start with an arbitrary element 



xo in X and define an iterative sequence by x n = T(x n -\) for n = 1, 2, 3, ... This sequence converges, 
and its limit is x*. The following inequality describes the speed of convergence: 



Equivalently, 



and 



d(x*,x n ) < - d(xi,x 

1 — Q 



d(x*,x n+1 ) < j—d(x n+1 ,x n ) 



d(x*,x n+ i) < qd(x*,x n ) . 



Any such value of q is called a Lipschitz constant for T, and the smallest one is sometimes called 
"the best Lipschitz constant" of T. 

Property 6. The fractional differential operators defined in (3 — 6) satisfy the following iden- 
tities: 

b b 

J f (x) Dl_g (x) dx = J g (xf £>«+/ (*) dx - f (x) I^g (x) \ b a 



or 





j f(x)D%_g(xf'D« + k(x) dx 

a 
b 

= J g (xf (xf D« + k (x) dx - f (x) I^g (xf D« + k (x) \ b a , 

a 

b b 

J f (x) D% + g (x) dx = J g (xf £>£_/ (x) dx + f (x) I^g (x) \ b a . 

a a 

Property 7. Assume a G (0, 1), j3 > a and / G C [a, b] . Then the following relations 



(8) 



(9) 
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= fix) 
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C Dl + I« + f(x) 


= fix) 


C D^_f(x) 


= fix) 



(10) 

hold for any x G [a,b]. Furthermore, the integral operators defined in (1,2) satisfy the following 
semi-group properties, 

jot r/3 _ T a+fS ja T P _ j-a+fS 

I a,+ I a,+ — V+ ' 1 b,- I b - ~ 1 b- 

Now, let's take up a fractional Sturm-Liouville problem for diffusion operator 



3. Main Results 

A Regular Fractional Sturm-Liouville Problem for Diffusion Operator 

Let's denote a fractional Sturm-Liouville problem for diffusion operator with the differential 
part containing the left- and right-sided derivatives. Let's use the form of the integration by parts 
formulas (8, 9) for this new approximation. Main properties of eigenfunctions and eigenvalues in 
the theory of classical Sturm-Liouville problems are related to the integration by parts formula for 
the first order derivatives. In the corresponding fractional version we note that both left and right 
derivatives appear and the essential pairs are the left Riemann-Liouville derivative with the right 
Caputo derivative and the right Riemann-Liouville derivative with the left Caputo one. 

Definition 8. Let a 6 (0, 1) . Fractional diffusion operator is written as 

L [Q;P;?] = D°_h (xf D% + + (2Xp (x) + q (x)) (11) 
where the function q (x) € L 1 [0, it] ,p(x) € L 2 [0, it] . Consider the fractional diffusion equation 

L ia,p,q]y>~ ( x ) + Xw <* ( x ) ( x ) = ( 12 ) 

where h(x) ,p (x) ^ 0, w a (x) > Vx S [0, n] and p, q, h, w a are real valued continuous functions in 
interval [0,7r]. The boundary conditions for the operator L are the following: 

ciyx (x) - c 2 ll;°h (Of Dl + y (x) = (13) 

x=0 

diy (tt) + d 2 ll;_ a h (nf D^ + y (vr) = (14) 

where d\ + d\ ^ 0, cf + c| ^ 0. The fractional boundary-value problem (12)-(14) is fractional 
Sturm-Liouville problem for diffusion operator FSLPDO. 

Teorem 9. Fractional diffusion operator is self-adjoint on [0, tt] . 

Proof: Let us consider the following equation 

TT TT 

{L[a, P ,q]<P,<f>) = J L[ a>p>q] <p(x) 4(x)dx = J <f>(x)D%_h(x) c Dg t+ <p(x)dx 



TT 

+ J (2Xp (x) + q(x))<p(x)<j)(x) dx. 



By means of property (8) and boundary conditions (13 — 14) and we obtain the identity 

IT 

(L[a, P , q] ^ <t>) = [ h ( x f D o,+^ (*f Do,+<P (*) dx + % (tt) cp (vr) + % (0) <j> (0) 



7T 

+ J (2Xp(x) + q(x))p(x)ej(x)dx. (15) 
o 

On the other hand, 

TT TT 

(f,L [a ,p, q }<P) = J P (x) L[a,p, q ]<P (%) dx = J if (x) D^_h (xf D^ + (p (x) dx 



TT 

J ip(x) (2Xp (x) + q (x)) (f) (x) dx. 



- 

'o 



Similarly, we obtain 

7T 

(<p, L [a ^ q]( f>) = [h (xf D% + <p (xf D% + cf> (x) dx + ^<p (tt) <j> (tt) + % (0) (0) 



7T 

+ j (2\p(x) + q(x))4>(x)ip(x)dx. (16) 


The right hand sides of the equations (15) and (16) are equal hence we may see that the left sides 
are equal that is 

{L[a,p,q\<P,4>) = {<P,L[a,p,q\4>) ■ 

Therefore, £[ a ,p,g] = L* apq y The proof is completed. 

Teorem 10. The eigenvalues of FSLPDO (12 - 14) are real. 

Proof: Let us observe that following relation results from Property (8) 

7T b 

j f (x) L [a , M g (x) dx = fh {xf Dg t+ f (xf D% >+ g (x) dx - f (x) l£°h (xf D^ + g (x) \% 

a 

7T 

+ j (2Xp (x) + q (x)) g (x) f (x) dx. (17) 


Suppose that A is the eigenvalue for (12 — 14) corresponding to eigenfunction y the following equal- 
ities are satisfy y and its complex conjugate y 

L [a, P ,q]y ( x ) + Xw °< ( x ) y( x ) = ° ( 18 ) 

ciyx (x) - c 2 ll~_ a h (Of D% + y (x) = (19) 
diy (vr) + d 2 ll;_ a h (irf D« + y (vr) = (20) 

L[ a ,p,q]y ( x ) + * w <* ( x ) y( x ) = ° ( 21 ) 

Cl y (x) - c 2 ll;_ a h (Of D% :+ y (x) =0 (22) 

diy (vr) + d 2 ll;_ a h (irf Dl + y (vr) = (23) 

where d\ + d 2 7^ 0,cf + c| 7^ 0. We multiply equation (18) by function y and (21) by function y 
respectively and subtract: 

(A - A) w a (x) y (x) y(x) = y (x) L [a , Ptq] y (x) - y (x) L [a , M] y (x) . 

Now, we integrate over interval [0, 7r] and applying relation (17) we note that the right-hand side 



of the integrated equality contains only boundary terms: 

7T 7V TV 

(A -A) J w a (x) y (x) y (x) dx = j y (x) L [a p q] y (x) dx - J y (x) L [apq] y (x) dx 



2\p (x) y (x) y (x) dx — J 2Xp (x) y (x) y (x) dx 
o o 

1-Q7 / \C 



l + y(x)I^h{xTDl + y{x) |J 



-y{x)ll : «h{x) C Dl + y{x) \% +y (x) I^h{xf D% + y (x) |J, 



-y{x)ll;«h{x) c Dl + y{x) 

IT 

(w a (x) + 2p(x))\y(x)\ 2 
o 

by virtue of the boundary conditions (19) , (20) , (22) , (23) , we find 

IT 

(A - A) J (w a (x) + 2p (x)) \y (x) | 2 dx = 
o 

and because y is a non-trivial solution and p (x) / 0, w a (x) > 0, it easily seen that A = A. This 
proves the theorem. 

Teorem 11. The eigenfunctions corresponding to distinct eigenvalues of FSLPDO (12 — 14) 
are orthogonal weight function w a on [0, tt] that is 

TT 

J (w a (x) + 2p (x)) y Xl (x) y X2 (x) dx = Ai ^ A 2 . 
o 

Proof. We have by assumptions fractional Sturm-Liouville for diffusion operator fulfilled by 
two different eigenvalues (Ai,A2) and the respective eigenfunctions {y\ 1 ,y\ 2 ) : 



L [a#,q]V\i ( x ) + Ai^ a (x) y Xl (x) = 



x=0 



ciy Xl (x) - c 2 ll^ a h(0) C D^ + y Xl (x) 
diy Xl (tt) + d 2 ll;_ a h (irf D% + y Xl (tt) = 

L [a,p,q]y\2 ( x ) + X 2W a (x) yx 2 (x) = 



x=0 



ciyx 2 (x) - c 2 ll;_ a h(0) C D« + y X2 (x 

diyx 2 (tt) + d 2 I 1 w ;_ a h (nf D* t+yx , (tt) = 
we multiply equation (24) by function y X2 and (27) by function y Xl respectively and subtract: 

(Ai - A 2 ) w a (x) y Xl y X2 = y^L^^y^ - y\ 2 L [a ^ q] y Xl . 



(24) 

(25) 

(26) 
(27) 

(28) 

(29) 



Integrating over interval [0, ir] and applying relation (17) we note that the right-hand side of the 
integrated equality contains only boundary terms: 

IT IT IT 

(Ai - A 2 ) J w a (x) y Xl (x) y X2 (x) dx = j y Xl L [a ^ q] y X2 (x) dx - J y X2 (s) L [a ^ q] y Xl (x) dx 
o oo 

IT IT 

= J 2X 2 p (x) y Xl (x) y X2 (x) dx - J 2A x p (x) y Xl (x) y : 



-y Xl (x)ll^h(xf Dl + y X2 (x) 

7T 

(Ai - A 2 ) J (w a (x) + 2p (x)) y Xl (x) y M (x) dx 



l + y X2 {x)I 1 -^h{xf Dl +VXl (x) |J 
-y Xl (x)ll;_ a h(xf Dl +VX2 (x) |5 
+y X2 (x)ll;_ a h(x) C D° + y Xl (x) ft. 



Using the boundary conditions (25) , (26) , (28) , (29) , we find that 

IT 

(Ai - A 2 ) J (w a (x) + 2p (x)) y Xl (x) y X2 (x) dx = 
o 

where Ai ^ A 2 . 

Let us now give certain auxiliary functions. Because we use the functions, the first of them is 
as follows 



\a-l 



(1,1) x - 

(a,-l) (a + 1,1) ~ vr-0 



where i^ 2 is the Fox- Wright function [13]. 

1*2 



(ai,ai) 

(Mi) (M2) 



£ 

fc=0 



r (ai + aik) z k 



the properties of the function are determined by the parameters 

Pi + @2 - «i = -1 



A 

S 

n 



l«i|" ai |/3il /3l l/3 2 l /32 = l 

1-2 „ 1 

&i + 6 2 - ai H — = 2a - -, 



(30) 



considering Theorem 4, we note that this function is continuous in [0,7r] when order a > 1/2, i.e. 
H > 1/2. For 0<a<l/2itis discontinuous at end x = it. The explicitly calculated function 
allows to estimate the second component of stationary function 4> of the differential part of diffusion 
operator 

D«_h(x) C D% +( f> (x) = 

which looks as follows 



^o (*) = £i + M + 



.7T — X) 



\a-l 



T(a)h{x) 



£ 1 +£ 2 il)(a,0,x) . 



(31) 



The next function is the following integral 

<p (x) = / V£-1 = ^+f^J = (* - 0) a (X - 0) a X! * 2 



(1,1) 

(a + 1,-1) (a + 1,1) 



x - 

7T - 

(32) 



Again, using the Theorem 4 and calculating parameters according to (30) 



we conclude 



A = -1, 5 = 1, [i = 2a + - 



a > => (i > 



and the obtained Fox- Wright function (32) is continuous in interval [0, ir] for any positive order a. 
Lemma 12. Let a > 1/2 and denote 

*A (y) = (2Ap (x) + g (x)) y x (x) + Aw a y A (x) 

A = cid 2 + c 2 cZi + c\diip (a, 0, 7r) 

Assume A / 0. Then, on the C [a, b] -space, the regular (12 — 14) is equivalent to the integral 
equation 

7T 

y x (x) = -/« + -J-_/«_Y A (y) + A(x)J Y x (y) dx + B (x) U + J^^-Yx (y)) 1*=* (33) 



where the coefficents A (x) and i? (x) are 



A (x) = [— c?2 + di^ (a, 0, x) — -0 (a, 0, 7r)] 
B (x) = ^ [ciV' (a, 0, x) + c 2 ] 

and functions tp is defined in (31) . 

proof. By aid of composition rules, equation (12) can be rewriten as follows: 



D«_h (xfDg, 



+ 



vx (x)+I^ + -^I^_Y x (y) 



0. 



The last equality refers that on the C [a, b] -space the function in brackets is a stationary function of 
fractional differential operator (FSLPDO). D" _h(x) C DQ + which according to (31) can be found 
as 

( \ce— 1 

<Po = £i + ^o Vr ( ~^ (x) = £i + ^ («, o, x) 

equation (12) in the form of 



yA (x) + J ° a+ h(x) ^- Yx {y) = il+ ^ (a ' °' x) ' 



(34) 



to end the proof we should connect coefficients j = 1,2 values c,-,dj j = 1,2 determining the 
boundary conditions (13 — 14). Let us note that the following formula results from composition 
rules (9) and equation (34) 

Il;_ a h (xf Dl +Vx (x) = -I^_Y X (y) + £ 2 , 



for continuous function y x , we obtain the following values a the ends 

Il;_ a h (xf Dq + i/x (x) = - [ Y x (y) + e 2 

x=0 J 



Il:_ a h(xf Dl +Vx (a 



e 2 



respectively for y x , using (34) we find 



y x (O) = O (0)=^ 



(35) 
(36) 

(37) 



y x (n) = fafr) - I« + —I*_Y x (y) 



= S 1 +Zrf( a ,0,„)-I° + —-I°_Y x (y) . (38) 

n \ X ) x=n 

The following set of linear equations for coefficients £j results from (35 — 39) 

c 2 X = c 2 £ 2 - (39) 
did + £2 {<h + diip (a, 0, tt)) = diF 

7T 

where ^ = J Y x (y) dx and F = I « + ^/£_Y A (y)\^. 


Since A^O, the solution for coefficients £j (j = 1, 2) is unique: 

Ci = {Fc 2 di - c 2 X (da + ciiV> (a, 0, tt))) /A 

£ 2 = (cidiF + c 2 d 1 X) /A 

substituting the above solution into (34) we recover the equivalent integral equation (33) . 
Furthermore we give notation such as 

A=\\A(x)\\, m p = min \p(x)\, 

xG[0,7rJ 

B=\\B(x)\\, M v = \\<p (x)\\ , 

where |.| denotes the supremum norm on space C [0, tt] . 

Teorem 13. Let a > 1/2. Suppose that A / 0. Then unique continuous eigenfunction y x for 
regular fractional Sturm-Liouville problem with diffusion operator (12 — 14) corresponding to each 
eigenvalue obeying 



\\(2Xp(x) + q(x)) + Xw a \\ < 



m r , 



M v + B<p (tt) + A (tt) m p 



(40) 



exists and such eigenvalue is simple. 

Proof. We have to say that equation (33) can be interpreted as a fixed point condition on 
function space C [0, tt] , 

Vx = Ty x 



where the mapping on the right-hand side is for any continuous function g G C [0, tt] defined as 

TT 

Tg = -^hjxj 1 ^ {v) + A ^J ^ (f ) dx + B ^ ( J o Q +7^^- y A (y)) U=, 



The following inequality will be useful in further estimations: 

\\Y x (g)-Y x (r)\\ = \\((2Xp (x) + q) g + Xw a g) - {{2Xp (x) + q) r + Xw a r)\\ 

= \\(2Xp{x) + q)(g-r)+Xw a {g-r)\\ = \\(g - r) ((2Ap (x) + q) + Xw a )\\ 
< \\g - r\\ \\(2Xp(x) + q) + Xw a \\ . 

By performing necassary operations for distance between images Tg and Tr for a pair of arbitrary 
continuous functions g,r € C [0, it] . 



\\Tg-Tr\ 



TT 

'^h^) 1 ^ +A ^j ^ (f ) ^ + ^ ^ {^h^) 1 ^ 



7T 

" /o " + 7^y (r) + A (x) / y * (r) dx + B {x) fa+hjz) 1 *'-^ (r) 



TT 

~ J ° Q+ (n (5) " Yx (r)) + ^ {x) J {Yx {9) " Yx (r)) dx + B (x) 



i 



-E_ (Y X (g) - Y X (r)) 



°' + h{ x y 

^ lb ~~ r ll • II (2Ap (x) + q) + Au> c 
< \\9-r\\L 



m. 



m r . 



where constant L = \\(2Xp(x) + q) + Aiy Q || + ^7r + ■ By means of (40) , we state that 

mapping T is a contraction on space (C [0, tt] , ||.||) 

||T 5 -Tr|| < || 5 -r||L L € (0, 1) . 

Hence, a unique fixed point enounced as y x G C [0, 7r] exists that solves equation (12, 33) and 
satisfies the boundary conditions (13, 14), provided (40) is applied. In that case such eigenvalues 
are simple. The proof is completed. 

Conclusion 

In this article, we have extended the scope of some properties of fractional Sturm-Liouville 
problem for the diffusion operator. We pointed that its eigenvalues related to the fractional Sturm 
Liouville problem for diffusion operator with the certain boundary conditions are real and its 
eigenfunctions corresponding to distinct eigenvalues are ortohogonal. Furthermore, we showed that 
the operator is self adjoint. These properties are similar to those for integer Sturm Liouville problem 
of diffusion operator. Our new results are original and important for the fractional Sturm-Liouville 
theory. 
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